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Abstract 

Electrons and electron neutrinos in the inner core of the core-collapse supernova 
are highly degenerate and therefore numerous during a few seconds of explosion. In 
contrast, leptons of other flavors are non-degenerate and therefore relatively scarce. 
This is due to lepton flavor conservation. If this conservation law is broken by 
some non-standard interactions, electron neutrinos are converted to muon and tau- 
neutrinos, and electrons - to muons. This affects the supernova dynamics and the 
supernova neutrino signal. We consider lepton flavor violating interactions mediated 
by scalar bileptons, i.e. heavy scalars with lepton number 2. It is shown that in case 
of TeV-mass bileptons the electron fermi gas is equilibrated with non-electron species 
inside the inner supernova core at a time-scale ~ (1 — 100) ms. In particular, a scalar 
triplet which generates neutrino masses through the see-saw type II mechanism is 
considered. It is found that supernova core is sensitive to yet unprobed values of 
masses and couplings of the triplet. 

Keywords: supernova, supernova core collapse, supernova neutrinos, lepton flavor vi- 
olation, bileptons, see-saw type II, scalar triplet, Higgs triplet model 



1 Introduction 



If the Standard Model (SM) is extended by a triplet of heavy scalar fields interacting 
with leptons and bosons of the SM, it becomes possible to generate Majorana neutrino 
masses through the see-saw type II mechanism [1]-|3]. Originally the mass of this triplet 
was assumed to be at a GUT scale; recently, however, TeV-scale see-saw became popular. 
In particular, it was argued that TeV-scale see-saw type II avoids electroweak hierarchy 
problem, encountered in GUT-scale see-saw, and meets the naturalness argument (see e.g. 
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jl][5]). See-saw II signatures at the LHC (see e.g. [6] [7] and references therein) and in 
low-energy experiments (see e.g. |S][H] and references therein) are extensively investigated. 

In the present paper we point out that TeV-scale see-saw II affects the core-collapse su- 
pernova physics, namely the supernova neutrino signal and the dynamics of the explosion. 
This conclusion is valid in a certain area of the parameter space of the neutrino mixing 
matrix. The effect is due to the lepton flavor violating (LFV) four-fermion interactions, 
which are induced by the scalar exchange. It is known that LFV processes may change 
drastically the conditions in the interiors of a core-collapse supernova [n][ID]lll Thus su- 
pernova physics appears to be an additional probe of TeV-scale see-saw II, along with 
LHC and low-energy experiments (such as rare lepton decays and muonium conversion). 

Scalar triplet, which enters the see-saw II, is a particular case of scalar bilepton, i.e. 
scalar which couples to a pair of leptons. Another type of bilepton is a doubly charged 
singlet scalar. If possessing TeV-scale mass, it may also be produced at LHC and influence 
the results of low-energy experiments. We include the analysis of its effect on the LFV in 
supernova in the present paper, because it mediates LFV interactions analogously to the 
doubly charged component of the scalar triplet. 

The paper is organized as follows. In section 2 we estimate the value of effective four- 
fermion LFV coupling necessary to significantly affect the supernova physics. In section 
3 we consider a simple extension of the SM with a doubly charged bilepton. Section 4 is 
devoted to see-saw type II extension of the SM. Influence of LFV inside the supernova 
core on supernova physics is discussed in section 5. The results are summarized in section 
6. 

2 Model independent considerations 

An explosion of a non-thermonuclear supernova is a manifestations of a collapse of a 
progenitor star core. Although numerous studies have failed to reproduce an explosion of 
a core-collapse supernova, there exists a commonly accepted general picture of the collapse 
(which is referred to as "standard picture" in what follows; for the details see, for example, 
|12j). When the mass of the iron core of a massive star reaches the Chandrasekhar 
limit, the infall phase of the collapse starts. The core contracts due to the gravitational 
attraction. Some fraction of electrons is converted to electron neutrinos through the 
inverse beta processes. When the density of the inner part of the core reaches the nuclear 
density value, ~ 3 ■ 10^^ g/cm^, the infalling matter of the outer core bounces from it. A 
shock wave emerges, which starts to propagate outwards. A powerful shock could expel 
the stellar envelope and produce the supernova explosion. However, according to the 
detailed numerical calculations [I3], the shock wave looses its energy and stalls after some 
tens of milliseconds. There exist several scenarios of presumably successful explosion (e.g. 
|16j). but neither of them is justified to date by detailed self-consistent calculations with 
neutrino transport. 

During the collapse extreme values of density, temperature, electron and electron 

""^In particular, in ref. |10| a Gelniini-Roncadelli model of neutrino mass generation [llj was proposed 
as a source of LFV. It resembles the see-saw type II model, but contains a massless Majoron, which is 
now excluded by the measurement of the Z-boson invisible decay width. 
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neutrino chemical potentials are reached. For the estimates we use the reference values 
of the supernova core parameters as presented in Table [H This data is obtained by 
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Table 1: Typical conditions in the inner supernova core (m < O.SM©) during the first 
50 ms after core bounce as obtained by us with the help of the open-code program 
Boom [T9j. Only SM interactions are taken into account. 

us with the help of the open-code programme [19] (in which only SM interactions are 
taken into account), being in good agreement with the results of more sophisticated 
calculations, see e.g. |2n]|2I]- The presented values are the mean values for the inner 
supernova core (m < O.SMq); the central values are even greater. Note that in the 
standard picture of the collapse electrons and electron neutrinos are plentiful and highly 
degenerate during the first hundreds of milliseconds. Their total number roughly equals 
the total number of electrons in the core before collapse, and decreases slowly due to the 
leakage of neutrinos from the core. By contrast, there are very few non-electron neutrinos 
and muons. Non-electron neutrinos have negligible chemical potentials; they are created 
in pairs with their antiparticles, concentrations being proportional to T^. Muons have 
chemical potential approximately equal to pe — Pu^ 5 their concentration has an additional 
suppression factor e~^™'^^~^'^^^'^ (see Appendix). Lepton flavor conservation law prevents 
electrons and electron neutrinos to turn into other types of leptons. 

This picture is modified if there exist LFV processes beyond the SM [5] [ID] • If electrons 
and electron neutrinos are converted into leptons of other flavors inside the inner cor^, 
the dynamics of the collapse and the supernova neutrino signal is different from the SM 
case. The following LFV reactions could be in principle relevant for supernovaH 



e e p p , 



(1) 



e^e^ — > eT p^ ■, 
e~Ve e"!/^,^, 

e'Ve P^Ve^r, (2) 

VeVe VeV^,^^, 
VeVe -> V^V^, 

The reactions are divided into two groups. In the first group all individual lepton flavor 
numbers, Lg, and are changed by two units or unchanged at all (e.g. for the 



^Not to be confused with ordinary resonant neutrino flavor conversion due to the MSW-effect, which 
take place in the envelope. It affects the neutrino signal but do not affect the supernova dynamics. 
•^We consider processes in which only leptons take part. 
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first reaction ALg = —2, AL^ = 2, AL^- = 0). In the second group at least some of 
tfie individual lepton flavor numbers are changed by one unit. As is discussed below, 
processes from the second group are tightly constrained by non-observation of rare decays 
and hardly can affect the supernova physics; in contrast, the processes from the first 
group can affect the supernova physics and at the same time fit the present experimental 
constraints. 

Let us estimate the value of the cross sections of the LFV reactions potentially impor- 
tant for supernova core. For this purpose we demand that such reactions should lead to 
the thermal equilibrium between e and v^. on the one hand and leptons of other flavor (-s) 
on the other hand at a certain timescale At (note that equilibration between charged lep- 
tons and neutrinos of the same flavor is established virtually instantly due to the ordinary 
weak interactions). We consider two reference values for At. The first one, Ati = 1 ms, is 
characteristic for the dynamics of the supernova (e.g. the shock propagates through the 
core during several ms). The second one, At2, should characterize the Kelvin- Helmholtz 
timescale of the cooling of a proto-neutron star (it is somewhat larger than the neutrino 
diffusion timescale, see section [5?T]) . which determines the long-term supernova neutrino 
signal. Although the inner core is deleptonized due to the diffusion at a timescale ~ 5 s 
|22j . we conservatively chose td^ = 300 ms, as at greater times the numbers presented in 
Table 1 are significantly changed. 

From the above-mentioned condition we find for the processes in which electrons are 
present in the initial state 

a > l/(Atnec) ~ 3 ■ 10-^^(300 ms/At)cm2. (3) 

If only neutrinos are involved in the process, the cross section should be several times 
greater, as the neutrino density is several times smaller. The estimated value appears to 
be rather small compared to ordinary weak cross sections. For example, 

a{uen e^p) ~ 2 ■ 10~^^ cm^ (4) 

for 160 MeV neutrinos [23]. 

Reactions ([T]) and ([2]) at low energies are described by effective four-fermion terms in 
the Lagrangian, the cross section being of order of G^jpynl^^, where He^ve. stands for the 
chemical potential of electron and electron neutrino correspondingly (every reaction is 
characterized, in general, by its own Glfv)- From ^ one gets a rough estimate for Glfv 
necessary for large LFV in supernova. In case of electron-electron collisions it reads 

Glfv > 4 ■ 10-V300 ms/At TeV~^ (5) 

In case of neutrino-electron and neutrino-neutrino collisions the bound is somewhat larger 
due to the smaller neutrino density and chemical potential. Such bounds on Glfv ~ A^/M^ 
indicate the mass M of intermediate heavy boson of order of several TeVs if the dimen- 
sionless coupling constant A is of order of 0.1. 

Are such values of Glfv compatible with present experimental constraints? This 
question is accurately treated in the following two sections, here we provide a preliminary 
discussion. The answer is different for two groups of reactions, ([T]) and ([2]). In the 
former case, the best experimental limit stems from the non-observation of oscillations of 
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Figure 1: ee — jjjj LFV transition mediated by the doubly charged bilepton A 

muonium atom, /i~e"^ <^ /^^e"^!]' reads Glfv ^ 10^^ TeV^^. Although muonium 
oscillations involve only charged leptons, due to the SU(2) symmetry the bound on the 
effective constant is valid for the first three processes in ([T]), including those which involve 
neutrinos. As for the fourth process in ([1]), it evades even this not very restrictive bound. 
Thus there is a room for LFV in supernova due to four reactions ([T]) with | ALg^^^rl = 0, 2. 

As for reactions with \ALi\ = 1 ([2]), they are tightly constrained by the rare lepton 
decays: 



The constraints span from Glfv < lO'^ TeV"^ (LFV decays of r) to Glfv < 10~^ TeV"^ 
— )■ e~e~e~^), see [25] and references therein. Therefore processes ([2]) can fit condition 
([3]) only marginally or can not fit it at all. Thus they are likely to be irrelevant for LFV 
in supernova. 

To conclude this section, we are interested in such an extension of the SM, in which 
four-fermion LFV processes with |AL;| = 1 are suppressed compared to the processes with 
|AL;| = 2. In the next section we give a simple but instructive example of such extension. 
In section 4 we show that see-saw type II model fits this requirements for certain values 
of neutrino masses and mixing parameters. 

3 Doubly charged scalar singlet 

Consider the SM extended by one doubly charged scalar A, singlet under SU(2). It can 
not be coupled to quarks due to its electric charge. The only possible coupling of such 
scalar to leptons reads 



Here 

/ and /' are charged leptons, /, /' = e, /i, r; 

A/// is a matrix of coupling constants, summation over lepton flavors is implied. 
C-conjugated right-handed fermions here and in what follows are defined as ip"^ = {iPrY = 
7^(V'k)*; standard representation of gamma- matrices being used. 

Let us choose the coupling matrix to be proportional to unit matrix: Xw = X6ii>. In 
this model low-energy processes with |AL;| = 1 are forbidden]^ 

^Note that in general the couphng matrix A = ||Ai;'|| is not invariant under rotations of the charged 



fi — )■ eee. 



(6) 




(7) 
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The exchange of doubly charged bilepton gives rise only to one LFV reaction, inter- 
esting in the context of supernova core conditions (see Fig. E]): 



e e ^ jjL jjL . (8) 
The corresponding effective low-energy LFV term reads 

^Iw = (9) 
which gives (see e.g. [26] [27] [28] ) 

cr(ee /i/i) = ^(1 - 2m^/s)|tmaa; - imm |/87r. (10) 

Note that one should take into account the combinatoric factors and the identity of 
fermions in the final state while calculating this cross section. In the center-of-mass 
reference frame \tmax ~ imin\ = ^E^-^ El — m"^. Taking Ef. c::^ fj,^ 2m^ one obtains 

.(ee^,ri-2-10-(^^j-^) cm^ (11) 
Comparing this with the estimate one finds that the supernova physics is modified if 

(MA/TeV)2 



> 10"V300 ms/At. (12) 



As was mentioned above, the best experimental bound on \^ / M\ follows from the non- 
observation of oscillations of muonium atom [2l]. We infer this bound from the review 
paper [5] taking into account that the definition of coupling constants in |5] and in the 
present work differ by the factor a/2: 

< 0.2, 90% CL. (13) 



(Ma/TcV) 



From inequalities f ll2p and (IT^ it follows that even for At = 1 ms there is a room for 
a singlet bilepton which affects the supernova physics and evades direct experimental 
constraints. If one takes A = 0.1, then the bilepton mass should satisfy 

220 GeV < Ma < 3.2 ■ (300 ms/At)"^/^ TeV. (14) 



lepton basis. It transforms according to A' ~ U'^AU ~ XU^U, where C/ is a matrix of a unitary rotation, 
WU = 1. However, in case when U does not contain complex phases one gets W = C/^ and A' = A. The 
role of complex phases is investigated in more detail in the next section. 
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4 Scalar triplet. See-saw type II. 



4.1 Cross sections 

Let us now discuss LFV mediated by a scalar triplet A. We are especially interested in 
the case when such triplet is responsible for the generation of Majorana neutrino masses 
through the see-saw type II mechanism [l]-[3]. In this mechanism scalar triplet addition- 
ally couples to the SM Higgs bosons, this coupling producing a vacuum expectation value 
for the neutral component of the triplet. The neutrino masses are proportional to this 
vev. 

The see-saw type II Lagrangian contains two major ingredients, a scalar-lepton inter- 
action, 

A/A = ^ \u'LfiT2ALif + h.c, (15) 
1,1' 

and a scalar potential, which in its minimal form reads 

V = -MlH^H + f{H^Hf + M^Tr(A"fA) + -^{fLH^ir2A^ H + h.c). (16) 

v2 



Here 



A° -A+/v^J' ^^^^ 



Li = ^ J is a doublet of left-handed leptons of flavor I = e, fi,T, H is a Higgs doublet, 

/i is a parameter with the dimension of mass. Note that the definition of coupling constant 
fl varies in different works by a numerical coefficient. Our definition coincides e.g. with 
that in [gj. 

Note that due to the anticommutation of the fermion fields 3x3 matrix A = ||A;;'|| is 
symmetric, 

= A. (18) 
The vev of the neutral component of the triplet reads 

(A°) = ^i-, (19) 

where v = \/2{H^) = 246 GeV. Due to the triplet vev neutrinos acquire the Majorana 
mass according to 

m = 2(A°)A, (20) 

where m = is a neutrino mass matrix in a fiavor basis. One gets that in the 

see-saw type II model neutrino mass matrix m is proportional to the coupling matrix 
A. Therefore now A can not be exactly diagonal as in the previous section, because m 
has non-zero non-diagonal entries (we know this since neutrinos oscillate). However, 
current experimental data on neutrino masses and mixing are incomplete and do not fix 
the ratio of diagonal and non-diagonal elements. In particular, it may appear that m is 
approximately diagonal, and thus it is possible to achieve the suppression of the processes 
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with |AL;| = 1 compared to the processes with jAL^I = 2. This issue is considered in 
detail in the next subsection. 

Lagrangian (IT3]) gives rise to all four reactions ([T]), potentially relevant for supernova. 
Their cross sections in the center of mass reference frame read 



a{ee ^ /i/i) = (|AeenA^^|7Mi)(l - m'j2Et)^l - my El E'j27r, (21) 

aieue ^ fxiy^) = i\\ee\^\\^^.\VMi){l - ml/AE^ fE^jATi, E, = E,^, (22) 
a{UeUe ^ U^U^) = {\Xee\^ \X ^,^,1^ M^) E^J 27r , (23) 

a{UeUe -> Url^r) = i\\ee\^\\rr\V M^) E^ /27T , (24) 

where neutrinos and electrons are assumed to be extremely relativistic. 

Density and chemical potential of electrons are higher than those for electron neutrinos; 
also an additional factor 1/2 persists in eq. f l22p . As a consequence, the first reaction in ([T]) 
dominates over the second and the third ones provided the c.o.m. energy is considerably 
greater than the threshold 2m^. The conditions on effective coupling constants, either of 
which ensures the appreciable rate of LFV in supernova core, now read 

AopA 



> 10-V300 ms/At (25) 
Ml ^ ^ ' TeV^ 



App A 



^ 2.4 ■ lO'VSOO ms/At ^. (26) 

It is shown in the following subsection that A^^ ^ A,-,- for the interesting set of neutrino 
parameters. Therefore it is always easier to satisfy condition fl2Sl) than condition fl2B]) : for 
this reason we concentrate on condition (125|) in what follows. 

If one assumes Aee — A^^ = 0.1 and takes At = 300 ms, then, analogously to the case 
considered in the previous section, supernova appears to be able to test Ma < 3.2 TeV. 
To compare, the LHC can probe Ma ^ 1 TeV [H]. 

4.2 Neutrino mass matrix, rare decays and experimental con- 
straints on TeV-scale see-saw type II 

LFV mediated by bileptons is constrained by non-observation of rare decays ([6]) and 
muonium conversion, see [5] [26]. The widths of rare decays and the rate of conversion 
are proportional to lAjj^jAigu/M^p, where l\ — I4 stand for the corresponding leptons. 
Experimental bounds are presented in Table [21 

Evidently the most stringent constraint stems from /i — )■ eee decay. In the present 
subsection we first concentrate on this constraint and address the question what mass- 
mixing pattern fits it and simultaneously allows sizable LFV in supernova. Then we show 
that if this constraint does not contradict the condition fl25l) . then other constraints also 
do not contradict it. 

The n — )■ eee constraint does not contradict the condition fl2^ if 

|A^^/Ae^| = |m^^/me^| > 40^300 ms/At. (27) 
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Table 2: Constraints from low energy experiments. We extract numerical values from 
the review paper [S] taking into account that the definitions of couphng constants in [5] 
and in the present work differ by the factor a/2. 

The neutrino mass matrix in the flavor basis, m, is obtained from the diagonal mass 
matrix through the transformation [29] 

m = U* ■ diag(mi, m2, ms) ■ U\ (28) 

where [/ = | |f/ij| | {I = e, fi,T, i = 1, 2, 3) is a PMNS neutrino mixing matrix: 

/I \ / ci3 srse-'' \ / Ci2 \ 

U = I C23 S23 \ [ 1 -Si2 Ci2 X 

V -S23 C23 / V -SlSe'^ Ci3 / V 1/ 

X diag(e'"^/2^gia2/2^^^)^ ^29) 

— COS ^ ij ; — SlXl ^ij • 

Note that eg. (1281) implies rri^ = m in accordance with eqs. (ll8l) . (l20|) . The explicit expres- 
sions for the entries of m may be found in [6l [71 [301 [31] j ^oi the sake of completeness we 
provide this expressions in our notationslfl 

ruee = « cfg + s23m3e2*'^ 

m^^ = mie-*"i (S12C23 + si3e-*'^ci2S23)^ + m2e-"^^ (C12C23 - si3e""^si2S23)^ + "^3C?3S23 
m^^ = mie-'"i (S12S23 - Si3e-*^Ci2C23)^ + m2e-*°2 (C12S23 + Si3e-*'^Si2C23)^ + m^cj^cl^ 
rrie^, = Cri[dsi2Ci2C23 + Si3e'^S23{m3- ae-'^'^)] 
rrier = Ci3[-rfsi2Ci2S23 + Si3e*''c23(m3 - ae-2*'^)] 

m^r = S23C23(-& + C?3m3) - Si3rfe-*''S12C12(C23 - S23) + Sf3ae-2"^S23C23 

(30) 

Here we define the parameters with the dimension of mass 

a = rriie'''^^ c\2 + m2e~^°'^ SI2 

b = mie-^"is?2 + m2e-'°2c?2 (31) 
d = m2e"*"^ — mie"*"^ 

^It is worth noting that if the Majorana phase matrix would stay at the first place in the r.h.s. of 
ea. ([^^ . then the absolute values of entries of m would not depend on Majorana phases. The correctness 
of the ordering of matrices in ea. (|^5)) is justified by the diagonalization of lepton mass matrices [25] . 
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If the masses mi and m2 are quasi-degenerate, mi ~ m2, and Majorana phases ai and a2 
are approximately equal, ai ~ 02 (this case is of particular interest, as is shown below), 
then this parameters reduce to 

a~6~e~*''imi (32) 

d ~ e-^"n^2 - mi) (33) 

There are nine parameters in the mass matrix: three neutrino masses m^, three mixing 
angles 6ij, the Dirac phase 6 and two Majorana phases ai^2- Experimentally measured (or 
constrained) quantities are [32i| : 

Ou ^ 34°, ^23 - 45" or 623 ~ 135", ^13 < 13°, (34) 

Aml^ ~ 0.8 • 10"^^ Am|2 ^ ±2.4 ■ lO'^eV^ (35) 
In addition, cosmological considerations [33] demonstrate that 

mi, m2, m3 < 0.35 eV. (36) 

Unknown quantities are phases 6 and ai^2, the mixing angle 6'i3, the sign of Am32, the 
quadrant of 6*23, the absolute scale of neutrino masses. Depending on the sign of Am^g two 
patterns of neutrino masses are possible: normal mass hierarchy (NH), mi < m2 < m3, 
and inverted mass hierarchy (IH), m3 < mi < m2. 

Formulas fl30p allow to check whether condition fl27p is valid for any specific choice 
of neutrino masses and mixing parameters. However, let us outline general trends of the 
dependence of |m^^/me^| on the neutrino masses, mixing angles and phases. 

Generically 

|m^^| ~ max{m2,m3}. (37) 

If one takes all three phases, ai,a2,S, randomly from the interval [0,27r], then in the 
majority of cases 

|me;,| ~ max{m2,si3m3}, (38) 

which implies |m^^/me^| < 5|^ Hence, in order to fit inequalities (l27|l one should first 
of all impose some bounds on the phases. Note that if all three phases vanish (which 
correspond to the exact CP-conservation in the neutrino sector), mg^ reduces to 

rrie^, ^ Si2Ci2C23Am^i/(mi + m2) + Si3S23Am^2/ ("^2 + "^3)- (39) 

The first term here is small compared to m2 if mi ~ m2. The second term is small 
compared to m3 if either m2 — m3 or S13 is small. Therefore we come to a conclusion that 
small phases, large masses (which implies degeneracy of masses) and small ^13 favour LFV 
in the supernova core. This is analogous to the conclusions of ref.[35], in which collider 
and low-energy LFV signatures were considered. 

We illustrate and refine this conclusions by means of plots depicted in Figs. [2ll6l On 
Fig. [2]the ratio |m^^/me^| is plotted as a function of a minimal neutrino mass (mi for NH, 

^To get this bound one should take into account that m2 > 0.009 eV and in case of normal hierarchy 
ms > 0.05 eV. 
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Figure 2: Ratio |m^^/me^| in case of exact CP-conservation (ai = 0^2 = 5 = 0). If 
^13 = then |m^^/me^| does not depend on the quadrant of ^23- Two thick horizontal 
hues correspond to \m^^/mf,^\ = 40 and |m^^/me^| = 720, which are minimal values 
sufficient for LFV in supernova core to proceed at scales 300 ms and 1 ms accordingly. 

ma for IH) for different values of ^13 and different choices of the quadrant of 623- All phases 
are taken to be zero, ai = 02 = S = 0. One can see that in case of quasi- degeneracy of all 
three neutrino masses, mi ~ m2 — m3 > 0.1 eV, condition fl27|) is fulfilled for any allowed 
value of 6^13. In case of a less stringent restriction, mi ~ m2 > 0.05 eV, condition fl27|) with 
At = 300 ms is fulfilled for 6'i3 < 3°. Note that two terms in rrie^ have different signs either 
in case of NH, C23 < or in case of IH, C23 > 0, see eq. (l39i) . This explains, particulary, 
the discontinuity at mi ~ 0.006 eV in case of normal hierarchy, ^13 = 3°, 623 = 135°. 

To what extent is it possible to relax the requirement ai = 02 = 5 = 0? To answer 
this question we first note the following facts. 

• If 

«! = 02, (40) 

then the parameter d reduces to e~*°^(m2 — mi), as was pointed out above. If we 
for a moment put ^13 = 0, then condition fHOl) ensures that me^ is generically small: 
rrief, = Ci3Ci2Si2C23e"'"i (m2 - mi) Am^i/ (mi + m2). 

• If a more restrictive condition is superimposed, 

«! = "2 = -2(5, (41) 

then rrie^ appears to be a sum of two generically small terms (with some relative 
phase factor) even for 6^13 7^ : 

rrie^, ~ Ci2Si2C23e^*'^Am^i/(mi + m2) + Si3e"^S23Am^2/ ("^2 + "^s)- (42) 

Condition (1411) is more general than the requirement ai = 02 = ^ = 0, and embraces 
the latter as a specific case. We illustrate allowed deviations from this condition on FiglHl 
One can see that equalities fHTl) should be valid up to ~ 3° for At = 300 ms in case 
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a2-«i 



Figure 3: Scatter plots for the ratio |m^^/me^|. Each point corresponds to a random 
choice of neutrino parameters in the following ranges: G [0.05 eV,0.35 eV],^i3 G 

[0°,3°],(5e [-90°,90°],ai G [-25-3°, -25 + 3°], ^2 G [-25-3°, -25 + 3°]. If the hierarchy 
is normal, then mmin = n^'i, otherwise mmin = ms 



when ^^13 < 3" and ^ 0.05 eV. One may further relax condition fj4T|) imposing more 
restrictive requirements on the degeneracy of neutrino masses and smallness of ^13. On 
the other hand the specific value of ai ~ 02 — —25, the quadrant of 6'23 and the hierarchy 
of masses generically do not affect the value of |wi^^/me^| significantly. 

Now it is easy to check that A^^ ~ A,-,- in the interesting range of neutrino mass- 
mixing parameters - a claim made in the end of the previous subsection. This statement 
is illustrated by FigHl 

Finally we should demonstrate that as far as condition (1271) is fulfilled, not only 
/X — )■ eee constraint but all other constraints from the rare decays are satisfied. Evi- 
dently, in order for the constraints cited in Table 2 to be compatible with condition (125 p 
one needs 

\mf,f,/mer\, \mf,f,/mf,r\, \mee/mer\, \mee/mf,r\, |"^eeMeAi| > 0.1. (43) 

Scatter plots for this ratios are presented on Fig|5l Evidently all the ratios under 
consideration generically satisfy the above constraint. However, the ratio |m^^/m^T-| 
appear to be smaller than 0.1 for a few choices of parameters. Let us consider this 
ratio separately in more detail. Taking into account that C23 = S23 one gets in the limit 
6^13 = 0, mi ~ m2 ~ m3 > 0.1, ai = 0^2 



+ m2e*' 



m3 — 17126 



(44) 



If «! = 02 = TT, this expression is small, \m^^/m^T-\ ~ |Am32/4m|| < 0.006. Therefore 
one should demand ai, 02 7^ tt in addition to eq. (HTj) in order not to fall in contradiction 
with r — )• /i/i/i experimental constraints. Numerical analysis (see FigJS]) shows that to 
evade the experimental constraint it is sufficient to have \ai — 7r| greater than few degrees. 
The above considerations demonstrate that the r — )■ /lee rare decay is the most useful 
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Figure 4: Scatter plots for the ratio Iw^/i/i/mT-rl- Each point corresponds to a random 
choice of neutrino parameters in the following ranges: Oi^ G [O^jS''],^ G [— 90°, 90°], ai G 
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one among all r — )■ III decays for the purpose of constraining the impact of LFV on the 
supernova physics. The importance of the r iiee measurement for probing the see-saw 
type II model was emphasized in [3l] . 



5 Possible consequences of LFV on supernova physics 
5.1 Standard collapse 

Before discussing the influence of LFV on supernova physics, we first give a sketchy 
qualitative picture using a standard approach (cf., e.g. j40J). 

There are several stages in standard core-collapse scenario where different neutrino 
processes dominate: 1) neutronization in the transparent core at infall before bounce; 
2) bounce, followed by shock-breakout with the most powerful neutrino burst; 3) accre- 
tion onto growing proto-neutron star with a stalled shock and possible shock revival; 4) 
Kelvin-Helmholtz stage of cooling of a hot proto-neutron star. The latter stage is char- 
acterized by a slow contraction of the hot core on a timescale of many seconds with the 
release of gravitational energy through neutrino emission. Contraction as a source of 
stellar luminosity was suggested by Kelvin and Helmholtz before the discovery of nuclear 
reactions. For the case of neutron cores inside stars this mechanism of stellar energy was 
considered by Landau [39] . 

For a polytropic equation of state 

P = Kp\ (45) 

the fundamental critical value of 7 is 4/3 in Newtonian theory. If the average 7 < 4/3 
over a stellar core, then this core is dynamically unstable. The only mass that may have 
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a hydrostatic equilibrium for 7 = 4/3 in Newtonian gravity is the Chandrasekhar mass 



This relation is derived in all classical textbooks on stellar structure. Physically, it is 
explained by a chain of simple arguments. The pressure P in general is the flux of 
momentum p, P pnv, where n is the number density and v is the velocity of particles. 
For extremely relativistic (ER) particles P ~ pnc, and for degenerate electrons p ~ 
PFermi ~ l/distauce ~ u^^^. By definition of Ye - the electron number density per baryon, 
n = pYe/rrin, where m„ is nucleon mass. This implies P oc n^^^ = {pYe/mnY^^ - which 
means 7 = 4/3 in relation fHS]) . (The definition of 7 is c^logP/Slogpl^.) 

We see that K oc {Ye/mnY^^ in (H5|) . On the other hand, gravity requires the equilib- 
rium pressure in the star center to be force/area oc GnM'^/R* oc GnM'^/'^ p^/^ . When 
equated with P = Kp'^/^ this fixes Mch in (H6l) . Since K oc {Ye/ninY^'^ one finds 
Mch oc (Fe/m„)^. When the coefficients are substituted accurately, one gets 



where mpi is Planck mass. 

In general, instead of Y^. there should be a sum over all particles producing pressure 
(e.g. e and v). And instead of the nucleon mass m„ there may be a mass of another 
particle if the latter dominates in producing gravity. E.g. if the Dark Matter would be 
dominated by a light fermion with mass mi the limiting mass ~ m^pjm\ would be on a 
scale of supercluster of galaxies if mi ~ 1 eV [50] . 

The contribution of baryons to the pressure is not significant until p > p„ ~ 3 x 10^^ 
(nuclear matter density). The equation of state during collapse is almost polytropic with 
7 slightly lower 4/3 due to small corrections from nuclei. The value of which depends 
on Yl = Ye + Yy and entropy S", is nearly constant since neutrinos are trapped. 

Of course, at the bounce and later, in the center of the core the pressure is dominated 
by nuclear forces. However, the initial collapse begins when the pressure is dominated by 
electrons. The onset of the collapse in standard picture is determined either by the drop 
of adiabatic exponent 7 below 4/3 due to iron dissociation or due to neutronization (or 
both, and in both cases the collapse begins due to the lowering pressure when the star 
core contracts). 

The dynamics of collapse in this case can be described by self-similar solutions [51] 
which show that a collapsing core separates into a homologous (f oc r ) central part, and 
an outer, supersonically infalling part. When the collapse is halted and bounces, a shock 
wave is produced at the sonic point. There is a self-similar solution with the shock as well 
[52]. 

A standard estimate for the free-fall collapse timescale [IQ] gives 




3/2 



(46) 





(47) 
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where pu is the average density in units of 10^^ g/cm^. Below we use notation for 
density in units of 10*^ g/cm^. 

This estimate is not very reahstic because it assumes zero pressure. In reahty the 
force produced by the pressure of ER degenerate leptons (~ n^^^R^) and the gravity force 
(~ GM'^ /E?) scale according to the same law (oc ). Therefore the contribution of the 
pressure force remains finite and non-negligible compared to the gravity force during the 
entire collapse. A better estimate for the collapse timescale tc which agrees with numerical 
simulations can be derived from self-similar models [52]: 

~ 6t/. (48) 

The mass of the inner core turns out to be about 10% greater than the Chandrasekhar 
mass for the same value of and is thus proportional to the lepton fraction Y^. 

Another important timescale characterizes neutrino diffusion. Only when the diffusion 
is less rapid than the collapse will the lepton fraction be frozen for the duration of the 
collapse. The standard expression for the diffusion timescale is 

td = ^^ . (49) 

crrin 

If one takes the values corresponding to the inner core already after bounce, M ~ 
O.6M0, p ~ pi4, a ~ 10"^^ cm^, one gets 

td ~ 10 s. (50) 

Kelvin-Helmholtz time is a bit larger than the diffusion timescale because must be 
multiplied by a factor reflecting the total energy storage over the energy in neutrinos. 

Neutrino signal basically consists of three consecutive periods. The first one is the 
neutronization burst - a short (~ 10 ms) powerful outburst of neutrinos which emerges 
when the shock wave reaches the neutrino sphere. Neutronization burst takes place few 
ms after bounce, and contains almost exclusively electron neutrinos, which are created 
during neutronization of stellar matter. The second period coincides with the existence 
of stalled shock and lasts ~ 1 s. The third period corresponds to the cooling of the proto- 
neutron star; it is characterized by an exponential decrease of neutrino luminosity and 
energy and lasts ~ 10 s. At the second and the third stages neutrinos and antineutrinos 
of all flavors are in approximate equilibrium at the neutrino sphere (not inside the core!) 
and therefore are emitted approximately in equal quantities. 

Neutrino oscillations take place above the neutrino sphere. In the layer ~ 100 km above 
the neutrino sphere collective oscillations occur (see e.g. [H]). In the outer stellar envelope 
more familiar MSW resonant flavor conversion take place (see e.g. [I5j). Oscillations 
modify the flavor content of the neutrino signal. In particular, they may convert a fraction 
of electron neutrinos of neutronization burst into non-electron neutrinos. 



5.2 LFV and collapse dynamics 

The rate of LFV reactions ([T]) during the infall phase is negligible, with possible exception 
of the last few milliseconds. There are two reasons for this. First, the cross-sections fl2Tl) - 
flMj) are small being proportional to the typical energies of e and Ue (i.e., to chemical 
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potentials of e and z/g). Second, the rates nf.^y^{^av) are small when the concentrations rze 
and TLy^ are low. Moreover, the first two reactions in ([T]) have thresholds proportional to 
muon mass m^. The reactions ([T]) may become important the few ms before bounce when 
the densities in the center approach the nuclear matter density. 

Let us estimate, for example, when ee — )■ /x/i process becomes important. The reaction 
is possible when > m^. At low density an electron may obtain such high energy only 
due to an exponential factor like exp(— m^/T), however T is always appreciably lower 
than m^. At higher densities electrons become degenerate, and we should estimate their 
chemical potential, taking into account that electrons are extremely relativistic: 

/ie^ll(pio>;)'/'MeV. (51) 

We see that /Xe surpasses the muon mass at p ~ 3 x 10^'^ g/cw? for = 0.3 and the 
process ee — )■ p/i sets in. However, if /Xg is only slightly above m^, only a small fraction of 
total electron density Ue = pYe/rrin contributes to nf,{av) of this reaction: only electrons 
near and above the Fermi surface take part in the reaction. 

Nevertheless, the process of transformation ee — )■ /i/i starts decreasing the electron 
pressure already at pis ~ 1, and becomes very important at pis ~ 10, when the chemical 
potential of electrons is twice the muon mass and the bulk of electrons contributes to the 
reaction. At this stage we may safely use an estimate for the reaction rate 

J ~ n,{av) ^ 300pi4ye s'S (52) 

where we take Ee = 2m^, A = 0.1, and Ma = 1 TeV (this gives o" ~ 2 x 10~^^cm^). Thus 
the timescale for the process ee — )■ pp is about 10 ms a few ms before the bounce. Our 
computations with the BOOM code show that the density p ~ 10^^ g/cm^ is achieved in 
the center about 5 to 7 ms before the bounce in good agreement with estimate (BSj) . and 
there is 1 to 2 ms from p ~ 10^^ g/cm^ until the bounce. 

The disappearance of electrons in ee — )• pp reaction starts actively contributing to the 
decrease of electron pressure - hence to the dynamics of the collapse. Though the effect is 
rather small (in a couple of ms we lose only 10 to 20% of electrons) still all small deviations 
of 7 from 4/3 may play a role. Moreover, the process of conversion will continue on later 
stages. It will lead to the decrease of electron (and electron neutrino) chemical potential 
and to the growth of entropy and temperature, since electrons are taken away well below 
the Fermi surface (cf. pTl ITS] ). 

It is interesting to estimate when the reverse process pp — ee becomes significant in 
order to speak about 7 in strict sense. See Appendix on equilibrium number of muon 
pairs in standard physics scenario when the muons would form a non-degenerate NR gas. 

The two-way LFV reaction ee ^ pp leads eventually to the establishment of equality 
p^ = pe for the chemical potentials. Muons become moderately degenerate since p^ — 
becomes somewhat larger than T. When the equilibrium is already established, the muon 
gas has 7 ~ 5/3 (P ~ pn^^v and v = p/m^, hence P ~ p^n^^/m^ oc n^J^ /m^)^ which 
is harder than for ER leptons, although softer than for nuclear matter of high density. 
However, in the beginning of the transformation ee — )• pp, when the reverse process is 
not yet effective, the pressure of the muon gas is much smaller than that of the electron 
gas due to the large mass (and therefore low momenta) of muons. Thus in the beginning 
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of the ee — )■ /i/x conversion, which may start several ms before the standard bounce, a 
decrease of Pg can not be compensated by the increase of P^, and we have an effect of 
reduction of total pressure similar to neutronization. This may lead to a faster collapse, 
hence to a stronger bounce. 

There is an opposite conclusion in [36]. They claim that the reduction of Yg (accom- 
panied by the reduction of Pg) leads to a weaker bounce. 

A one-zone collapse model was used in [35]. It may be good for a uniform homologously 
contracting core where 7 is close to 4/3. Even when electrons are transformed into 
trapped electron neutrinos the gamma for lepton gas is close to 4/3 (all leptons are 
extremely relativistic). In this framework they give a standard estimate of the mass of 
the homologous core (HC), which is just approximated by a Chandrasekhar mass Mch '■ 
Mhc — Mch — 5.8Y^Mq (a bit more accurate estimate following [SH SD] is Mhc — 
l.lMch). 

However, if the core started collapsing and then drops fast enough and non- 
uniformly (e.g., faster in the center, as in our case), then the picture may be more 
complicated. 

Effective 7 lower than 4/3 means faster collapse, hence a partial loss of homology. 
Quantitatively, the effect will be small if the process ee — )■ fifi is very slow. Anyway, here 
we see the qualitative difference with the standard case which may lead to the deformation 
of the homologous flow at infall and even to the formation of a second strong shock 
(because the innermost parts of the flow, those that are in the deepest potential well, have 
to accelerate due to the decrease of pressure). We postpone the discussion of quantitative 
results for future research keeping in mind possible new dynamical phenomena. 

Another place where the reduction of 7 may be interesting is the shock after the 
bounce. There is a stage of more or less steady accretion through a stalled shock. When 
the pressure in the shock is dominated by leptons, and not by nuclear matter, the reduction 
of 7 means higher compression (the density jump in any strong shock is (7 + l)/(7 — 1). 
However, when the free path corresponding to LEV processes ([T]) is much larger than a 
hundred of km this effect will not be noticeable. 

Above we mostly concentrated on the LEV conversion between electrons and muons. 
A related questions should be answered if electron neutrinos equilibrate with muon and/or 
tau neutrinos. How does this affect supernova dynamics? Does this help the supernova 
to explode, or prevent the successful explosion? 

A comprehensive analysis of similar problems may be found in paper [3S]. A lepton 
number violating reaction UeN — )■ UeN is considered there. As in [37] it is found that the 
decrease of lepton degeneracy leads to a weaker shock wave, which prevents the prompt 
explosion (which is anyway disfavored by numerical simulations). Another consequence 
is the increase of entropy and temperature of the core. In principle this could affect the 
delayed explosion scenario. However, the authors of [SB] state on the basis of numerical 
simulation that the effect of LEV on the shock wave revival is negligible. A difference 
between [38] and present work should be pointed out: Ue -H- z/g transition results in a 
zero chemical potential of electron neutrinos with no compensation in other neutrino 
flavors, while LEV transitions considered here produce degenerate seas of muon and/or 
tau neutrinos at cost of reducing chemical potentials and concentrations of electrons and 
electron neutrinos. Eor this reasons the results of [5B] can not be directly applied for LEV 
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processes. 



5.3 LFV and supernova neutrino signal 

First let us consider whether modifications of the neutronization burst are possible. Neu- 
tronization burst consists of neutrinos sweepped out by the shock wave from the outer 
core. On the other hand, as discussed above, the LFV reactions ([T]) may proceed effectively 
only in the inner core where the density and lepton chemical potentials are sufficiently 
high. Therefore it does not seem plausible that the fiavor content of the standard neu- 
tronization burst may be altered. However, as pointed out in the previous subsection, 
LFV may modify the shock wave itself, which may cause the modification of the form, 
duration and total energy of the neutronization burst. If additional shock waves are pro- 
duced due to inhomologous collapse, they result in additional neutronization bursts, in 
striking contrast to the standard picture. 

Now let us turn to the neutrino emission after the neutronization burst. Its spectrum 
and fiavor content is determined solely by the equilibrium conditions on the neutrino 
sphere. How can the redistribution of chemical potentials between different flavors inside 
the inner core affect such emission? To answer this question numerical simulations are 
necessary; here we propose only one possible signature. LFV leads to lowering of the 
average neutrino energy; also a fraction of electron neutrinos (which take part both in 
neutral current and in charge current interactions) is converted to non-electron neutrinos 
(which take part only in neutral current interactions). Both this circumstances tend to 
increase mean free path of neutrinos and to accelerate the cooling of the core. On the other 
hand, a decrease of neutrino degeneracy tends to decrease mean free path and to decelerate 
the cooling. An interplay of this opposite trends may modify the total duration of the 
neutrino signal. This signature may be easily revealed by modern neutrino detectors, 
which should detect thousands of neutrino events in case of galactic supernova explosion. 

6 Summary 

We have shown that the lepton flavor violation, generated by the exchange of TeV-scale 
scalar bileptons, may drastically modify conditions inside a core-collapse supernova. An 
important example of such bileptons is a scalar triplet responsible for the neutrino mass 
generation in the see-saw type II scenario. In this scenario the matrix of scalar-lepton cou- 
plings is proportional to the neutrino mass matrix. The latter is currently partially flxed 
by neutrino oscillations and cosmological data. Also the effective four-fermion couplings, 
generated by the triplet, are restricted by low-energy experiments. We analyze the body 
of data and demonstrate that there exists an allowed region of the see-saw II parameter 
space in which LFV is strong enough to lead to the thermal equilibrium between elec- 
tron and non-electron species inside the collapsing core. Roughly speaking, LFV affects 
the supernova physics if the matrix of scalar-lepton couplings is approximately propor- 
tional to the unit matrix, Xui ^ XSw, and the effective coupling A^/M^ is greater than 
10~^ TeV~^ ~ 10~^Gf- The former condition corresponds to neutrino masses greater than 
0.05 eV and neutrino phases which satisfy the relations \ai + 26\, \a2 + 26\ <^ n. 
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Equilibration between electron and non-electron species is a striking feature, absent 
in the standard picture of the collapse. It starts to settle not earlier than few milliseconds 
before bounce, when the density exceeds 10^^ g/cm^. Such equilibration lifts the chemical 
potentials of muons and non-electron neutrinos higher than 100 MeV, somewhat decreases 
the chemical potentials of electrons and electron neutrinos, leads to the increase of tem- 
perature and entropy. In addition to Fermi seas of extremely degenerate electrons and 
electron neutrinos, which are present in the standard picture, seas of extremely degenerate 
non-electron neutrinos and moderately degenerate muons emerge in case of LFV. This 
may affect the supernova physics in various ways. An elaborated study should be carried 
out in order to accurately and thoroughly describe all possible consequences. In particu- 
lar, collapse simulations with new interactions included are necessary. In the present work 
we restrict ourselves by just outlining various possibilities. They include modifications of 
supernova dynamics and modifications of supernova neutrino signal. 

Possible modifications of supernova dynamics are basically due to the decrease of pres- 
sure in the center of the core, where LFV reactions proceed effectively. Such decrease gives 
rise to a non-homologous stage of the collapse which may start few ms before the "stan- 
dard" bounce. This may lead to modification of the shock wave and even to additional 
shock waves. In addition, the increase of temperature and entropy may affect the delayed 
explosion scenario. 

Possible modifications of supernova neutrino signal include: 

• Modification of the form, duration and total energy of the neutronization burst, 
which accompanies the shock breakout. If additional shock waves emerge due to 
LFV, then additional neutronization bursts are expected. 

• Modification of the total duration of the supernova neutrino signal, which signifies 
the alteration of the cooling rate of the proto-neutron star. There are three effects 
which affect the neutrino diffusion from the inner core: the emergence of seas of 
degenerate non-electron neutrinos, which do not take part in charge current inter- 
actions; the decrease of the average neutrino energy and the decrease of neutrino 
degeneracy. However, the first two effects tend to speed up the cooling, while the 
third one - to slow down. In the present work we do not find which trend wins. 

The discussed signatures in neutrino signal do not concern the flavor content of the 
neutrino signal and therefore are not smeared by neutrino oscillations. 

We emphasize once more time that all effects listed above should be verified by pro- 
found supernova calculations. Once such calculations are done, it is possible to use the 
results in various ways. On the one hand, one can try to construct a successful supernova 
explosion model taking into account lepton flavor violation. On the other hand, one can 
extract parameters of the see-saw type II model and other models which include bileptons 
(or restrict such parameters) from the already persistent or future astrophysical data. In 
particular, one may get some bounds from the very fact that supernova explodes. 

As a finale remark we note that the discussed signatures of LFV in supernova are rather 
general and may be relevant in the context of various microscopic models, not necessarily 
those which include bileptons. In particular, the required LFV may be generated by 
the spin flips of Majorana neutrinos with large magnetic moment [H] or by the flavor 
changing neutral currents in SUSY with R-parity violation [42] . 
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Appendix: Comments on number of muon pairs in 
collapsing supernovae 

The role of muons in collapsing supernovae was pointed out long ago by Domogatsky 
[13] . He has shown that for T > 20 MeV the number density of muons is sufficient to 
trap muon neutrinos u^^ due to charge current reactions with muons only. The role of 
neutral currents is understood much better nowadays, and we know that u^^ and must 
be trapped even without muons at lower temperatures due to scattering on nuclei and 
electrons. The role of muons is important in the models considered in this paper, and 
also, e.g., for leptonic photons discussed by Okun [H] (see also |15], [16] and references 
therein) . 

It is trivial to find the number of muon pairs for the case of standard physics following, 
e.g.. Landau & Lifshits [47j). Here we add just a few methodological remarks (and point 
out one place where Landau & Lifshits make an incorrect statement). 

The photons inside stars have the blackbody Bose distribution, their chemical potential 
being equal to zero. The muons and antimuons are in thermal equilibrium with the 
photons, particularly, due to the reaction 



H + photons , 



(53) 



and this means that [17] 

where /i^- and are the chemical potentials for /i~ and respectively. 
On the other hand, reactions 



(54) 



e +p 



z/e + n 



(55) 



ensure that 



^^^^ f^i^ij. ~ A'-e /^i/e- 

It may be shown that in the supernova core fii, is negligible, which leads to 



(56) 



1^)1 — f^e l^Ue- 



(57) 
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Boom code simulations show that /i^ does not exceed 60 MeV. 
The number of muons is 



r 47rp^ _dp_ 

go = 2 being the statistical weight of a muon. We are interested in the conditions when the 
temperature is of order of a few MeV or of a few tens MeV, so muons are non-relativistic 
and non-degenerate, then one can use the Boltzmann-like approximation: 



goeM-im, - ^.,)/T) ^^^^^y^^^j.^ " ^ = ^^o exp(-(m,-;.,)/T) 




Let us put yU^j = for a moment, then an exponential factor reads exp( 
At the first glance it seems surprising that the number of pairs is expressed through 
exp(— m^/T) and not through exp(— 2m^/T), while the energy required to create a pair 
is 2m n- A simple explanation of this paradox can be suggested using the Boltzmann 
distribution. Let us denote the number of /x"*" in the volume V as A^"*", then indeed, 
relative to vacuum {Nq = 1) 

^ = ^exp(-2m,/T) (60) 

where statistical weight g~^ must take into account the contribution of all translational 
states of the muon /x"*" in the volume it occupies V/N~^ = l/n^- The integral over 
d^p/ {2u:Y with the Boltzmann factor exp(— p^/2m^T) gives (m^T/27r)^/^. The same factor 
is supplied to g^ by each born together with the /i"*" somewhere in the volume V , so 
this must be multiplied by the whole V now, and we find 




(61) 

Since A^"*" = A^~ we recover n^n~ found above in (1591) . 

It is interesting to compare the value of from fl59p with the electron concentration 
rif,. We must notice that in the book by Landau & Lifshits |17| §105 it is written not 
correctly that for T > me the number of electron-positron pairs is high. It is true only 
for the "terrestrial" densities. Inside the stars the situations when T > me, but T < fie, 
are not rare. The equilibrium condition for electrons does not imply the equality of the 
electron and positron concentrations: when the temperature is lower than the chemical 
potential He, the number of the positrons is exponentially lower than that of the electrons 
by a factor of exp(— 2/ie/T). During the collapse, the value of density is like p = 2 x 10^^ 
g/cm^ and for Ye = 0.3 this gives fj,e = 200 MeV in eq. (l5Ti) . so the number of positrons 
is small even for T of order tens MeV. The correct expression for electrons is 

ne = £-,f^l. (62) 
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Thus in the standard scenario 



V2 



exp(-(m^ -/x^)/T) 




(63) 



in the supernova core. We see that LFV dramatically raises the number density of muons 
making ~ ng. 
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